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Abstract 

The theoretical foundations of a new general approach to the measurement 
problem of vibrational observables in trapped ion systems is reported. The 
method rests upon the introduction of a simple vibronic coupling structure 
appropriately conceived to link the internal ionic state measurement outcomes 
to the mean value of a motional variable of interest. Some applications are 
provided and discussed in detail, bringing to light the feasibility and the wide 
potentiality of the proposal. 
PACS: 03.65.Ta; 32.80.Pj; 42.50.Ct; 42.50.Hz 
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I. INTRODUCTION 



The question of how measuring the probability distribution of the eigenvalues of an 
observable A in a given quantum state often exhibits, at least, a high level of practical 
complexity. Sometimes, however, all what we want to know about the physical quantity 
under scrutiny is just its expectation value. To extract this information is of course possible 
if the probability for each eigenvalue of A in the given state is known, even if, in principle, 
such a detailed knowledge goes beyond the specific need. 

In the last few years advancements in ionic trapping and cooling techniques have provided 
the possibility of carefully controlling the dynamics of a two-level object coupled to a bosonic 
system 0,0] ■ I n a Paul microtrap, for example, a suitable configuration of inhomogeneous 
and alternating electromagnetic fields force the centre of mass of a confined two-level ion to 
behave as a quantum harmonic oscillator @,|]. When appropriate laser beams act upon the 
system, a vibronic coupling is switched on. Thus one may efficiently manipulate the trapped 
ion quantum state in order to realise, for example, applications in the quantum logic context 
and non-classical states 0. 

Proposals of measurement processes aimed at revealing the quantum dynamics of such 
systems have been reported. 

In 1997 L.G.Lutterbach and L.Davidovich have presented a proposal to extract all pos- 
sible information about the vibrational quantum state of a trapped ion centre of mass 
Such a knowledge is reached via the reconstruction of the Wigner distribution, W(a), rela- 
tive to the ion centre of mass degrees of freedom (a = x + ip, x and p being the dimensionless 
harmonic oscillator position and momentum respectively). It is well known that the value 
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of the Wigner function in the phase space origin coincides with the mean value of the parity 
operator in the relative quantum state M. Moreover, subjecting the system to some suitable 
external forces whose action may be described by displacements operators, one may induce 
a phase space translation of the physical system and hence a phase space translation of the 
Wigner distribution, passing from W(ce) to W(a) : W(a — a ) = W(a), a being related to 
the external action. Exploiting the two features just recalled one has, in principle, a recipe 
to measure W(cx) at any point. Indeed the mean value of the parity operator evaluated 
in the system Wigner function translated by a^, is just the value of the original Wigner 
function in the point a^. 

If the interest is pointed in the individuation of the vibrational quantum energy prob- 
ability distribution only, one can turn to quantum nondemolition (QND) measurements as 
already proposed in two-level trapped ion systems both in one- dimensional P,|10|| and two- 



dimensional [II] contexts. Generally speaking, in QND measurements a signal observable 
is indirectly measured observing a meter quantity which has been suitably modified via 
an appropriate signal-meter interaction [12,13|. If the signal is conserved both by the free 



hamiltonian and by the signal-meter interaction, then the signal observable may be repeat- 
edly measured getting rid of disturbances due to the "free time evolution" (between two 
measurements) or to the "back-action" (during the signal- meter interaction). 

In the trapped ion context, QND schemes have been applied taking the vibrational en- 
ergy as the signal and the atomic populations as the meter [|9|-|TT|j. The method is based 



upon the ability of implementing a suitable vibrational energy conserving vibronic coupling 
that resonantly drives weak atomic transitions of the atom beyond the Lamb-Dicke limit 
(signal-meter interaction), and upon the possibility of monitoring such transitions (meter 



observations). For example in ref [10 it has been proposed a quantum nondemolition mea- 
surement scheme of the vibrational population of a trapped ion, by means of Raman pulses 
inducing resonant electronic transitions. The fundamental point of this method is the de- 
pendence of the Rabi oscillations between the two electronic levels on the energy quantum 
number. 

The two protocols just recalled allow to realise measurements of quantum observables 
or even reconstructions of the entire dynamical configuration of the system, but both are 
generally very complicated to be executed. On the other hand, in many cases the knowledge 
of the mean value of some suitable observable is enough to reveal the phenomena we are 
interested in. 

In this paper we present a new general approach to the measurement of the mean value 
of a vibrational dynamical variable in the trapped ion context. The method, its reliability 
and its feasibility are described in detail in the second section. In the subsequent section 
some possible applications are described bringing to light the effectiveness and the wide 
potentiality of our proposal. 



II. THE IDEA AND THE METHOD 

The physical system on which we focus consists in a two-level ion confined into a three- 
dimensional Paul trap |L]-|4j]. In such a trap the motion of the ion centre of mass is confined 
in an effective harmonic potential which, under appropriate conditions, may be rendered 
isotropic |Q, so that the free dynamics of our system is governed by the free hamiltonian 

Ho = huiT X! + h^A&z (1) 
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Here lot is the frequency of the trap, Oj(a|) with i — x,y,z are the annihilation (creation) 
operators of the centre of mass oscillatory motion, ua is the atomic transition frequency and 
a z the z component of the pseudospin operator a= (a x ,a y ,a z ) associated to the internal 
ionic dynamics. Here a x = |+) (— 1 + |— ) (+|, o y = (— | — |— ) (+|), a z = |+) (+| — |— ) (— | 
and |±) are the two effective atomic levels such that a z |±) = ± |±). 

In the experiments currently performed on trapped ion systems, the probability of find- 
ing the ion in its ground (|— )) or excited (|+)) electronic state may be measured using 
quantum jumps techniques with almost 100% efficiency JL4]. The measurement of bosonic 



dynamical variables of interest is in practice systematically traced back to atomic population 
observations. 

In this paper we sketch an experimental scheme that, in principle, allows to measure the 
mean value of vibrational dynamical variables. The key idea of our proposal relies on the 
ability of implementing an appropriate vibronic coupling by which we succeed in making 
the atomic populations after the interaction dependent on the mean value of the vibrational 
observable in the initial state, that is before the interaction. In other words, starting from a 
dynamical configuration, we induce a time evolution able to correlate electronic populations 
to the mean value of the observable in the initial vibrational state. 

Let A be the dimensionless vibrational (or bosonic) observable under scrutiny, given 
in the Schrodinger picture, and assume it is time independent. Suppose now to subject 
the system to a suitable laser configuration responsible for realising a vibronic coupling of 
strength 7 that, in the interaction picture relative to Ho, has the form 

HP = h 1 Aa x (2) 
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We will show that the hamiltonian given in eq.(@) is responsible for inducing atomic 
transitions modulated by the mean value of A, < A >, evaluated in the initial state. 
In the simple case of a trapped ion prepared in the state 

10(0)) = \a) I-} (3) 

with \a) such that A\a) = a \ a), the interaction picture quantum dynamics under the action 
of the hamiltonian given in eq.(|]) is 

i/)W(tf) = cos(7«t) | a) |— ) - isin^yat) \a) |+) (4) 

Looking at this time evolution we note that the Rabi oscillation frequency is proportional 
to a. Hence, as given in advance, the flipping rate from |— ) to |+) is just proportional to 
the mean value of A in the initial vibrational state: ("0(0)1 A |0(O)) = (ipmbrl A \ip V ib r ) = ex. 

Obviously in the more general case of an arbitrary bosonic initial state 

1^(0)) = Nw) I-) (5) 

the analysis is not so simple. Nevertheless it is always possible to establish a correlation 
between electronic transitions and the initial mean value of A. To this purpose, subject 
the ion to a |-pulse relative to the interaction of the form h5a x (6 being a real number 
representing the strength of the interaction), easy to implement with a laser tuned on the 
electronic transition frequency (ua) HJH- ^ ne effect on the system is the passage to the 
state 

m) = ^ mb r)\-)y (6) 

where |— ) is an eigenstate of the operator a y and |VWw) is the same generic vibrational 
state appearing in eq. (^j) . 
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Calculate now the time evolution of the a z mean value under the action of the hamil- 
tonian in eq.(Q). Since [H ,a z } = the interaction picture operator coincides with the 
Schrodinger picture operator a z , and hence 



(*,(*)) = ^ [P [I) it)) = (^\t)\a z \^\t) 

'm) u\(t)d z u A {t) m 



(7) 



with 



u A {t) 



-i-4-t 



(8) 



the interaction picture time evolution operator. 
Moreover, being 



(j\{t) a z U A (t) = cos(2-yAt)a z + sm(2<yAt)a y 



(9) 



one immediately obtains, for the state in eq.(O), 



(a z (t)) = - (if) vibr \ sm(2jAt) \ip vibr ) 



(10) 



For the sake of completeness, consider the case of an operator having a non-degenerate 
spectrum being partly continuous and partly discrete. Let VL D = {a n } be the set of all the 
discrete eigenvalues and flc the set corresponding to the continuous spectrum of A. 

The vibrational state \tpmbr) can be expanded in terms of the states \a) as follows: 



\lpvibr) = ( a \ ^vibr) |«) da + 



aeflr 



aeQc 



(a\ ip V ibr) | a) da 



(11) 



Starting from this expansion, the mean value of sin(27At) appearing in eq. (|T0|) , may be 
cast in the form 
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(Vwl sin (2 7 At) \i> vibr ) = \( a \ i'vibr) | 2 sin (27 at) da + 

/ I (oi\ i^vibr) | 2 ^ai{2^at)da (12) 

In many physical situations of interest, it is reasonable to assume that | (a\ i/) V ib r ) | 2 may 
be neglected when a is outside the interval [— a max , ct max \, a max being an appropriate positive 
real number (we shall refer to this assumption as finite effective spectrum hypothesis) . At 
this point we wish to emphasise that the possibility of such an assumption is not necessarily 
related to an a priori knowledge of \ip V ib r )- Consider for example the case of A = X (X 
being the x-component of the ion centre of mass position operator given in the Schrodinger 
picture). For such an operator it results Qd = and Qc = R- Nevertheless, due to the ion 
confinement, one may be confident in the fact that all the effectively involved values of the 
observable X must be compatible with the trap dimensions. 

As another example consider the bosonic excitation number operator hi (with i — x, y, z). 
One has f2c = and, in principle, £l D = N. Nevertheless, again because of the confine- 
ment, the probability of finding the charged trapped particle with rii beyond an appropriate 
maximum value, may be reasonably assumed as being negligible. 

It is worth noting that the z component of the orbital angular momentum L z , is canoni- 
cally equivalent to the operator h x — h y . Thus, bearing in mind the previous considerations 
regarding hi, we legitimately claim that also L z satisfies the finite spectrum hypothesis in 
the experimental situation under scrutiny. We remark that this property holds for any 
component of the orbital angular momentum due to the spherical symmetry of the prob- 
lem. Analogous considerations are applicable to the correlation operators as, for example, 
C xy = (a^cby + dyCix), reaching the same conclusions. Indeed, like L z , also the operator C xy 
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is canonically equivalent to n x — n y , being e % * Lz C xy e~ % % Lz = e l ^ Cxy L z e~ l ^ Cxy = h x — h y . 

Under the finite spectrum hypothesis we may control both the coupling parameter 7 and 
the interaction duration t in such a way that sin(27«t) ~ 27«t for any a inside [— a max , a max \. 
(In the following we shall refer to the interval of real values for which sinx ~ x as the 
linearis ability zone.) 

Such an assumption makes it legitimate linearising the mean value of the operator 
sin(27yit), which means, reducing it to the mean value of 2^ At. Indeed, since the only 
relevant contributions to the mean value calculations are the ones relative to a belonging to 

Q! c = Q c n [-a max , a max ] or Q' D = Q D n [-a max , a max ], one deduces 



(ifj vibr \ sin(2~f At) \ifj vibr ) = I ( a \ ^vibr) | 2 sin(27at)da 

+ I (a I ip V ibr) | 2 sin(2 / yat)da ~ 
— X! \ (d\ ip V ibr) \ 2 sm(2 j at) da + 

aeQ' D 

+ I (a\ ip V ibr) | 2 sin(2 / yat)da ~ (13) 

Since the calculation is developed in the linearizability zone, it is possible to linearise the 
sinusoidal function before coming back to the original dominion of the variable a: U 
Hence it follows that 

{ipvibr\ sm(2^fAt) ~ I (a\ ifj vibr ) \ 2 2~fatda + 

aeQ' D 

+ / I (a I Vw) \ 2 2~fatda ~ 
~ Y I (a| Vw) \ 2 2~fatda + 



+ I (a I ip V ibr) \ 2 2'jatda = 

Jaefln 



(ip vibr \2-fAt\ip vibr ) (14) 



After such a discussion we are ready to deduce the main result of this paper. Assuming 
that A has an appropriately bounded spectrum and accordingly controlling the pulse area 
jt, we may take advantage of eq. (|I4D into eq.(|IO"|) obtaining 

(A(t = 0)) = (ij mbr \ A \ifj mbr ) ~ -— (a z (t)) (15) 

Eq. fll5D relates in a transparent way the mean value of the vibrational observable A 
before the interaction in eq.(^) (that is in the initial vibrational state), with the mean value 
of the z-pseudospin component, (a z (t)), measured immediately after this interaction. 

We again emphasise that the applicability of eq.(|15"D requires, among other conditions, 
the narrowness of the interval containing the eigenvalues of A appearing in the expansion 
of the initial vibrational condition. The importance of this point may be fully appreciated 
considering that the identification of the mean value of sm(2jAt) with the mean value of 
2jAt amounts at neglecting the mean values of all the nonlinear terms appearing in the 
series defining the operator sin(2j At), that is at assuming 

Indeed, the finite spectrum hypothesis allows to claim that | (A n ^ \ < a^ax legitimating 
relation ( fL6|) when ^a max t « 1. 

For the sake of clarity we summarise the steps of the process we propose for measuring 
the mean value of vibrational observables in the state \ip V ibr)'- 

a) Estimate an upper limit, a max , for the modulus of the bosonic eigenvalues of A involved 
in iV'wbr)- As previously stressed, such an estimation does not require the a priori knowledge 

Of \tpvibr}- 
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b) Convert the initial condition given by eq.(5) into the state \ip) = |Vw>r) \ ~) y by the 
action of a |-pulse relative to the interaction h8a x . 

c) Subject the system to a pulse related to the interaction = ^A&x. The intensity and 
the duration of the pulse must be chosen in such a way that 2 r )a max t lies in the linearizability 
zone. 

d) Measure the atomic populations, calculate (a z ) and hence deduce {ip V ib r \ A \ip V ib r ) from 
eq-flUD. 



We conclude this section giving some orders of magnitude for the parameters involved 
in the measurement process. To this end we consider for example A — hi. According to ref. 
|IJI4|, Hi = 20 is surely a good estimation for a maximum value of the bosonic excitations, 



in the sense that the experiments up to now performed generally involve just few phonons. 

About the coupling strength 7, we may refer both to the Innsbruck group experiments 
(with Ca + ions), where typical 7 values are of the order of magnitude of lOKHz [TE\, or 



to Boulder group experiments (with Be + ions), where typical 7 values are of hundreds of 



KHz |]14 |. Although the latter values are higher than the former ones, we guess it should 
be possible to obtain smaller values even with Be + ion experiments just decreasing the laser 
intensities. 

Assuming therefore 7 ~ lOKHz and [—0.4,0.4] as the linearizability zone, in the case of 
A — fii, one estimates an interaction time t ~ lfj,s which is in the grasp of experimentalists 
PJ. Indeed we emphasise that the brief duration of the vibronic interaction is compatible 
with the laser "turning-on"-"turning-off"and also makes our method reasonably immune 
from decoherence effects. 

11 



III. APPLICATIONS AND CONCLUSIVE REMARKS 



A fundamental underlying requirement of the method presented in the previous section 
is the implementation of the hamiltonian given by eq.(0) relative to the specific bosonic 
observable one wishes to measure. Fortunately the trapped ion context allows to tailor 
practically at will vibronic couplings possessing a desired form. By changing indeed the 
lasers configuration and parameters, a huge class of effective Hamiltonian operators may be 
engineered. 

In this section we focus our attention on some vibrational observable variable of physical 
interest bringing to light the effectiveness and the wide potentiality of our method. In 
particular we shall concentrate our attention on those vibrational operators of clear physical 
meaning, for which the possibility of applying the hypothesis of finite spectrum has already 
been discussed in the previous section. Taking in mind the conclusions reached in section II, 
we now therefore limit ourselves to show the feasibility of the coupling mechanisms required 
in order to apply our measurement procedure. 

a) First of all let's consider the bosonic excitation number fv In order to measure the 
mean value of such an operator exploiting the method presented in this paper, the effective 
Hamiltonian model 

H\ I] = h-yfii&x (17) 



has to be implemented. In ref [16] it is shown how this model can be realised in laboratory. 



b) Another important quantity to measure is the angular momentum of the centre of 



mass of the confined particle. In ref. |17]], in connection with other purposes, it has been 
implemented the interaction hamiltonian 
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H { / ] = hjLi&x (18) 

where Li is the component of the orbital angular momentum along the direction i = x,y, z. 
Such an interaction is useful to reach our goal with respect to the angular momentum 
operator. 

We remark that, at the best of our knowledge, our proposal provides for the first time a 
method exploitable for measuring this observable in trapped ion systems. 

c) The next operator which we wish to apply our measurement protocol to, is a space- 
correlation operator defined as 

Cij = Sja] + hjCil = QiQj + PiPj i ^ j (19) 

with i, j — x,y,z and where Qt = ^(^fc + %) an< ^ ^ k = ~ ^k) are the dimensionless 

ion center of mass position and momentum operators respectively. 

The interaction hamiltonian model of the form (0) corresponding to A = may be 
easily implemented. Consider for example two Raman schemes [I] responsible for two cou- 
plings realised via two effective lasers directed along the axis x' and y', rotated by f with 
respect to x and y axis. Suppose that the two effective lasers relative to each Raman scheme 
are tuned on the carrier frequency u>a, are n out of phase and finally have wave vectors of 
equal moduli. Under these hypotheses the total interaction hamiltonian takes the form 

a x (20) 

where a± = ^= (a x ± a y ), r\ is the Lamb-Dicke parameter and / an assigned function giving 
rise to non-linearity effects stemming from the non vanishing gradient of the electric field 

ma- 
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in 



f(r],a 1 + a + ) - ffadLh-) 



If we restrict our analysis to the Lamb-Dicke limit (that is assuming n << 1), up to 
terms of order r/ 2 , we just obtain 

Hp = hjC xy a x (21) 

with 7 oc Trj 2 . 

Measuring the operator C xy is interesting, because, for example, it provides a simple way 
for revealing the so called "Parity effect" predicted the first time in the context of bimodal 
cavity QED [18] and then in the context of two-dimensional trapped ions In ref. [19 



it has been proved that a trapped ion, prepared in a SU{2) coherent state and subjected to 
the action of a suitable set of lasers whose effects are described by a bimodal and two-boson 
Jaynes-Cummings model, H\ = hj (d x a y |+) (— | + h.c), exhibits a nonclassical behaviour 
that may be brought to light measuring the mean value of C xy . This operator, indeed, 
distinguishes in a macroscopic way between two just microscopically distinguishable initial 
states. 

More in detail, it has been analytically predicted that, if Nq is the total number of bosonic 
excitations in the initial state then, at particular time instants t* and t*, it results 

, l+A^o if N odd 

C xy ) (t* p ) = (22) 

| —Nq if N even 

+N 2 if N odd 



respectively. 



(23) 

— Nq if N even 



The possibility of measuring the correlation C xy and the operator C 2 y is then useful to 
observe the "Parity" discrimination. 
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d) As last example we consider the dimensionless position operator Q x = -h^(a x (X - 
X. Implementing the hamiltonian 

H { P = 7 (a x + al)a x (24) 

is very simple. Indeed it requires the use of two effective lasers both directed along x and 
suitably tuned. One of the lasers must be first red sideband tuned and the other one must 
be tuned on the first blue sideband. With such a laser configuration one obtains 

= 7l (a, |+) (-1 + at |-) (+|) + 72(4 |+) (-| + a x |-> (+|) (25) 

Adjusting the lasers intensities to have 71 = 72 = 7, the interaction hamiltonian model 
coincides with the one expressed by eq . ([24]) . 

Summarising, in this paper we have discussed a simple method aimed at measuring the 
mean value of vibrational observables. The key idea of our proposal is the possibility of 
transferring information about the expectation value of a vibrational operator to the atomic 
populations, with the help of a suitable vibronic coupling. Applications discussed in this 
section highlight the relevance of the present procedure to a wide variety of situations. 
Apart from some specific applications here considered, we emphasise that our method is 
also compatible with non-spherical traps. The reason is that our proposal rests upon the 
implementation of the interaction in eq.@, which does not depend on the trap isotropy. 

Concluding we wish again to stress that, when we want to know only the expectation 
value of a physical quantity in an assigned vibrational state, the method presented in this 
paper might be easier to realise with respect to the others already quoted in the introduction. 
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